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Recent experimental progress in coupling nanoscale conductors to superconducting microwave 
cavities has opened up for transport investigations of the deep quantum limit of light-matter inter- 
actions, with tunneling electrons strongly coupled to individual cavity photons. We have investigated 
theoretically the most basic cavity-conductor system with strong, single photon induced non-local 
transport effects; two spatially separated double quantum dots (DQD:s) resonantly coupled to the 
fundamental cavity mode. The system, described by a generalized Tavis-Cummings model, is inves- 
tigated within a quantum master equation formalism, allowing us to account for both the electronic 
transport properties through the DQD:s as well as the coherent, non-equilibrium cavity photon state. 
We find sizeable non-locally induced current and current cross-correlations mediated by individual 
photons. From a full statistical description of the electron transport we further reveal a dynamical 
channel blockade in one DQD lifted by photon emission due to tunneling through the other DQD. 
Moreover, large entanglement between the orbital states of electrons in the two DQD:s is found for 
small DQD-lead temperatures. 

PACS numbers: xx 



I. INTRODUCTION 

In circuit quantum electrodynamics (QED), the meso- 
scopic analog of cavity QED, solid state qubits are cou- 
pled to superconducting microwave cavities on chip 
Circuit QED systems combine the appealing properties 
of high cavity quality factors and strong vacuum mi- 
crovifave fields with low qubit decoherence. This has al- 
lowed for experiments in the strong coupling limit with 
qubit-cavity coupling exceeding the qubits decoherence 
rates. The strong cavity-qubit coupling together with 
fast, coherent manipulation of the qubits has lead to an 
astonishing development in the areas of quantum infor- 
mation processing 0-0] and microwave quantum optics 
with superconducting circuits (8l-[l3j. Moreover, circuit 
QED architectures have a large potential for simulations 
of strongly interacting many-body system s llill a nd tests 
of fundamental quantum physical effects |15l4l7| . 

The rapid development in circuit QED triggered in- 
vestigations on nanoscale qubits and conductors cou- 
pled to microwave cavities or resonators [nil. Par- 
ticularly interesting are recent experiments on few-level 
quantum dots coherently coupled to microwave cavities 
IM m m [13-14^. These experiments open up for 
transport investigations of light-matter interactions in 
the deep quantum regime: single electrons interacting 
strongly and coherently with individual microwave pho- 
tons. The large versatility of microwave photon state 
properties [l^ |4l|, , together with the well-established 
controllability of quantum dot levels provide a broad 
scope for fundamentally important experiments. 

A key feature of conductor-cavity systems is the possi- 
bility to coherently couple electrons in conductors sepa- 
rated up to centimeters [40[ . This puts in prospect entan- 
gling macroscopically separated transport electrons, of 
importance for nanoscale quantum information process- 
ing and Bell inequality tests fist - lioj . Moreover, this non- 
local feature can be harnessed for efficient heat-transfer 



or refrigeration over large distances [47|, |48| . A first step 
towards these goals would be an experimental demonstra- 
tion of non-local, few-photon mediated, electronic trans- 
port effects. In the present work we investigate theo- 
retically the simplest possible strongly coupled cavity- 
quantum dot system where such non-local effects can be 
observed: two double quantum dots (DQD:s) coupled to 
the same transmission line cavity (See Fig. [T]). We argue 
that dot-cavity systems with single-level or metallic dots 
will, in comparison, display strongly surpressed non-local 
effects. 




FIG. 1: Two DQD-cavity system: Each DQD (pair of red 
ovals) is coupled to the central conductor of the transmission 
line cavity (middle blue rectangle), two lead electrodes (gold 
rectangles) and to two gate electrode potentials (silver rect- 
angles). Here V^L(H)i and Vgi^(^R)i denote the left (right) lead 
and gate electrode potential of DQDi, respectively. 



The DQD:s are resonantly coupled to the same, fun- 
damental mode of the microwave cavity. This DQD- 
cavity system constitutes a generalized Tavis-Cummings 
model [4^, [H^l, with strong hybridization of the DQD 
electron and microwave photon states. Our investigation 
is focused on the non-local electronic transport prop- 
erties. For a broad range of parameters, the current- 
voltage characteristics provide clear signatures of trans- 
port electrons exchanging photons. Further, for asym- 
metric DQD-lead couplings, noise and higher order cur- 



2 



rent fluctuations reveal a dynamical channel blockade in 
one DQD lifted by single photons emitted by electron- 
tunneling in the second DQD. In addition, we demon- 
strate the existence of large orbital entanglement between 
electrons in different DQD:s. 

We emphasize that the predicted non-local effects are 
direct consequences of the non-equilibrated, transport- 
induced photon state. This makes our investigation qual- 
itatively different from earlier transport studies on pairs 
of two-level systems coupled via thermalized bosons [sil . 
[5^ . To fully account for the coherent, non-equilibrium 
properties of the photon state, as well as the electron 
tunneling through the DQD:s, our investigation is carried 
out within the framework of a quantum master equation 
(QME). The approach is similar to the ones used to inves- 
tigate transport through single two-level swtems coupled 
resonantly to a photonic mode in Refs. psi. [23. Is^]. 

The article is organized as follows: In Sec. Ill Al we 
introduce the model for the closed DQD-cavity system, 
present its Hamiltonian and discuss the eigenstates and 
eigenenergies. We further derive, in Sec. Ill Bl a QME 
for the system when the DQD:s are coupled to lead elec- 
trodes. In Sec. mil the transport properties are investi- 
gated in the regime where the temperature of the lead 
electrodes exceeds the DQD-cavity coupling strength. 
Focus is put on the non-local transport properties, cal- 
culating the non-local current- voltage characteristics, the 
current cross-correlations and the full counting statistics. 
In Sec. II VI we turn to the regime where the lead temper- 
ature is smaller than the DQD-cavity coupling strength. 
Transport signatures of coherent electron-photon inter- 
action as well as entanglement between the electrons in 
the DQD:s are investigated. The effect of dephasing and 
approaches to minimize this effect are discussed in Sec. 
lYl 



II. SYSTEM AND METHOD 

The system considered is depicted in Fig. [T] Two 
DQD:s, denoted 1 and 2, arc inserted near the cndpoints 
of a transmission line cavity. The central conductor is 
capacitively coupled to the right (left) dot in DQD1(2) 
(see e.g. Ref. [37| for a possible experimental realiza- 
tion) . One gate and one lead electrode is further coupled 
to each dot in the DQD:s. The leads are assumed to be 
in thermal equilibrium with a common temperature T 
and chemical potentials ^^i, with v = L,R and i = 1,2 
denoting to which dot the lead is coupled. 

Throughout most of the paper we will consider the 
strong coupling limit. This implies that the DQD-cavity 
coupling is large compared to the DQD-lead couplings 
and also dominates over decoherence due to other type 
of system-environment interaction. Moreover, we will the 
DQD-lead couplings to be much stronger than the inter- 
action with the rest of the system-environment and hence 
neglect decoherence from the latter. Only in last section 
the effect of DQD dephasing as well as DQD-rclaxation 



and photon loss will be considered. 



A. Model 

The DQD:s, forming singly occupied two-level systems, 
couple linearly to the microwave photons in the cav- 
ity. The system Hamiltonian Hs, describing the DQD- 
photon interaction as well as the orbital degrees freedom 
of the DQD:s and the direct interaction between the DQD 
charges, is derived in Appendix A. Below we will take the 
DQD:s to be on resonance with the fundamental mode of 
the transmission line cavity. Moreover, the cavity char- 
acteristic impedance Zq is assumed to be much smaller 
than the resistance quantum Rq = h/e^, relevant for 
regular transmissionline cavities. Under these conditions 
the DQD-cavity system will be described by a generalized 
Tavis-Cummings (TC) Hamiltonian [4^ 
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Hs = hojo'a + 

i 

+ hgo{d'' dl^dei + h.c.) 



(1) 



Here and {dgi and dei) denotes the creation (an- 
nihilation) operators of the ground and exited, i.e. of 
the bonding and anti-bonding, state of DQDi. We have 
further introduced the photon creation operator and 
the frequency lj of the fundamental mode, and the DQD- 
cavity coupling strength go, for simplicity taken equal for 
both DQD:s. Note that since we have assumed single oc- 
cupancy of the DQD:s, the spin-degree of freedom of the 
DQD:s will only have the effect of renormalizing tunnel- 
ing rates and is hence neglected in Eq. ^ and below. 

The generalized TC Hamiltonian in Eq. ([1]) has the 
form of a TC Hamiltonian for both DQD:s occupied while 
it reduces to a Jaynes-Cummings (JC) and a harmonic 
oscillator (HO) Hamiltonian when one or none of the 
DQD:s are occupied, respectively. It is of key impor- 
tance for the discussion below to describe the eigenstates 
in the HO, JC and TC subspaces of Hs- We first note 
that Hs commutes with the operator for the number of 
excitations n = a^d+X^i ^L^ei- The eigenstates can then 
be characterized by the corresponding quantum number 
n. We express the eigenstates in terms of the DQD-cavity 
product states \£,i£,2p), with DQDi in the state with 

= 0, g, e, and p photons in the cavity mode. 

For the HO subspace £,1 — £,2 = and hence the num- 
ber of excitations is equal to the number of photons, 
giving the eigenstates |00n). In the JC subspace with 
DQDI (2) occupied the eigenstate with zero excitations 
is \S^0) = \gOO) i\S^O) = \OgO)). For states with a fi- 
nite number of excitations the photon state and the state 
of the occupied DQD hyridises. The states, denoted by 
\Sfn), are superpositions of product states with n and 
71 — 1 photons in the mode. For DQDI occupied they 
arc given by \S^n) — [\gOn) ± \eOn— l)]/^/2 and for 
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DQD2 occupied we have {S^n) = [|05n)±|0en- l)]/\/2. 
The eigenbasis in the TC subspace has a similar struc- 
ture. The state with zero excitations is a product state 
|D°0) = \ggQ) and the states with one or more excitations 
are superpositions of product states with different num- 
ber of photons. We denote the finite-excitation eigen- 
states by \D°'l), with a = 0, ±, for n = 1 and \D'^''n), 
with and /3, 7 — ±, for n > 2 and give their exact forms in 
Appendix B. The spectra of the HO, JC and TC Hamil- 
tonians, also given in Appendix B, are shown in Fig. [51 
Importantly we see in Fig. [5] that a state with n excita- 
tions has an energy nhuj + O{hgo) relative to the energy 
of the state with zero excitations in its subspace. More- 
over the TC (HO) ground state is shifted ~hw/2 {fujj/2) 
with respect to the JC ground states. 
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FIG. 2: The spectrum of the Hamiltonian in Eq. IH) with the 
eigenenergies marked by the horizontal black lines. Column to 
the left shows the spectrum of the HO subspace (both DQD:s 
unoccupied), middle column shows JC spectrum (DQDl or 
DQD2 occupied) and right column shows the spectrum of 
the TC subspace (both DQD:s occupied). The diagonal red 
arrows show the allowed transitions due to electron tunneling 
described by Eq. ((3|. 



B. Quantum Master Equation 



The leads and their tunnel coupling to the DQD:s are 
described by the Hamiltonians and Ht, respectively. 
The lead Hamiltonian reads Hi = ^ ^ ^kC^^^Cki^i, with 

c^^j creating an electron in the state with energy in 
the lead connected to dot vi. In the ground-excited basis 



the tunneling Hamiltonian is given by 



kLi 



s\\i{9i)dgi + cos{6i)de 



tmcl 



km 



cos{9i)dgi + sm{9i)de 



+ h.c], 



(2) 



with t^i denoting the energy independent lead-dot tun- 
neling amplitude for dot vi. We have further introduced 
the DQD mixing angles tan(0i) = ±tifli/(-\/w2 — ± 
w), where denotes the interdot tunneling amplitude 
of DQDi. The -l-(-) sign here refers to the energy differ- 
ence between orbitals of the left and right dots of DQDi 
being positive (negative). 

We assume weak tunnel couplings between the dots 
and the leads and restrict the investigation to the se- 
quential tunneling regime. Following the standard Born- 
Markov approximation scheme a quantum master equa- 
tion (QME) is derived for the time evolution of the re- 
duced density matrix p of the DQD-cavity system jssj . 
We point out that the lead-dot tunneling rates T^i = 
X]fe "^(^ ~ ^fc) must be chosen much smaller the 
DQD-cavity coupling strength go- This restriction is nec- 
essary for the strong coupling condition to hold. More- 
over, it allows us to neglect coherences between states 
with an energy difference Ae > hgo, i.e. to perform a sec- 
ular approximation. Considering for simplicity identical 
tunnel-couplings to the left and right dot in each DQD, 
i.e. = Tin = r,, we can write the QME ^ = £[p], 
with the Liouvillian 



C[p] = --[Hs,p] 



EE 



d€de'r^^i{6i) 



U^{e)d^^5{e + e' - Hs)dU{e - Hs)p 



+Lrie)dyie + e' + Hs)d^^6ie' + Hs)p 
-fAe)5{e + e' + Hs)dy{e' + Hs)pk^ 

"/.,:(e)'5(e + e' - Hs)d.,S{e' - Hs)pdl + h.c. 



(3) 



Here f^i{e) = /(e - p^i) and f^i{e) ^ 1 - /(e - p^i), 
with / denoting the Fermi-function. We have further 
introduced the rates TLeiidi) = ^Rgi{Si) ~ TiCos^{9i) 
and TlM = fRe^{e^) = r,sm^i9,). 

We point out that electron tunneling into or out of the 
DQD:s in the sequential secular regime can be associated 
with jumps between energy levels in adjacent columns in 
Fig. [2l This means that the most compelling physi- 
cal picture will emerge when working in the eigenbasis 
of Hs (see Sec. Ill A[) . By evaluating the rates for the 
different tunneling processes, described by Eq. in 
this eigenbasis we can draw the following two important 
conclusions: 

(i) The number of excitations of the system can change 
with or — (-|-)1 when an electron tunnels into (out of) 
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one of the DQD:s, as indicated by the red arrows in Fig. 
[5J The system energy can thus increase or decrease by 
the energy fujj/2 + O{hgo) in a tunnchng event. Conse- 
quently, an electron tunneling through one of the DQD:s, 
e.g. from the left lead into the DQD and then out to 
the right lead, can change the system energy by 0, ±ficj 
+O{hgo). If a tunneling electron changes the energy with 
+(—)huj + O{hgo) a subsequently tunneling electron, in 
the same or the other DQD, can absorb (emit) this en- 
ergy. This process, here referred to as the transport elec- 
trons exchanging a photon, is depicted in Fig. |3l 

{ii) The rate for processes where the energy is in- 
creased or decreased for an electron tunneling into dot 
v in DQDi is proportional to T^eif^i[huj/2 + O{hgo)] 
and T^gif^i[~huj/2 + 0{hgQ)], respectively. Similarly, the 
rates for exciting or deexciting the system when electrons 
tunnel out of DQDi are given by r^gif^i[—haj/2+0{hgo)] 
and T^eifmi—hhj /2 + 0{hgQ)]. Note that the rates for 
tunneling events in which energy is emitted or absorbed 
can be controlled by the chemical potentials and mix- 
ing angles di. Note further, that for go = these rates 
coincide with the corresponding rates for the DQD:s de- 
coupled from the cavity mode, as depicted in Fig. |31 

Based on an analytical solution to the QME (See 
Appendix C), we can further conclude that a well de- 
fined steady-state solution for p exist for mixing angles 
6*1, 02 > 7r/4. For smaller angles and certain bias config- 
urations the photon number can diverge. In Ref. [3l|, 
where a single DQD-cavity system was considered, it was 
shown that DQD mixing angles 6 < tt/A can lead to pop- 
ulation inversion and hence a cavity lasing state. Since 
the focus on this work is the few-photon regime we, if 
not otherwise explicitly stated, focus on mixing angles 
01, 02 > 7r/4 below. 



III. NON-LOCAL TRANSPORT PROPERTIES 

The main purpose of this work is to investigate non- 
local transport properties due to exchange of photons 
between tunneling electrons in different DQD:s. This in- 
vestigation is carried out in the regime huj k^ T ^ %o- 
The experiments reported in Refs. [H, [H, were per- 
formed under these conditions. Moreover, for ksT ^ hgo 
the occupations of the leads do not change significantly 
over the energy scale hgg, i.e. f^i[zthaj/2 + 0{Hgo)] w 
f,^i{zthuj/2). Then there are effectively only two ener- 
gies in the leads at which the electrons can tunnel into 
and out of the DQD:s. This simplifies the expressions 
for the tunneling rates into and out of the DQD:s (See 
Sec. Ill Bp and allows to reduce the QME to an ordinary 
master equation (ME) (See Appendix C). We write this 
ME dP/dt = MP, where P is a vector with the prob- 
abilities of eigenstates of Hs and M is the matrix with 
the transition rates between these eigenstates. 

For the system to display non-local transport effects it 
is clear that two conditions must be fulfilled: First, trans- 
port electrons in different DQD:s must exchange photons. 
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FIG. 3: (a) Elementary process in which two transport elec- 
trons exchanges a photon. In step I an electron tunnels 
through DQDI, i.e. into (blue arrow) and then out of (red ar- 
row) DQDI, while increasing the system energy by hoj. In step 
II an electron tunnels through DQD2 and changes the system 
energy by — ftoj. Through steps I and II the two electrons have 
exchanged a photon (whiggling line) through the cavity, (b) 
Tunneling processes into and out the ground and excited state 
of DQDi decoupled from the cavity mode. The rates for the 
processes 1,2,3, and 4 are f Le.ifLi{fii^ f LgifLi{—huj/2), 
fReifRi{Fi'^/2) and f RgifRi{—hw/2). (c) Bias configuration 
of DQD2 considered in Sec. IIIII The positions of the chem- 
ical potentials of the left leads /iL2 and /i_R2 relative to the 
energies at which electron tunnel into and out of the DQD 
are shown. These sets of energies are in an interval ~ hgo 
denoted by the grey boxes around the dashed lines. 



Second, for at least one of the DQD:s the effective tun- 
neling rate into the empty DQD or out of the occupied 
DQD must be dependent on the number of excitations 
in the system or the occupation of the other DQD. The 
most clear non-local effects will thus occur for a state- 
dependence such that transport becomes blocked in one 
of the DQD:s if photons are not emitted by the other. 
Here we take the DQD2 to be blocked when there are 
no excitations in the system. This is accomplished by 
choosing a bias configuration similar to the one depicted 
in Fig. O i.e. such that /l2(-^w/2) = 1, fR2{fuj/2) = 0, 

fL2{-huj/2) = /L2(fi^/2) = 1. 



A. Current 

We first consider the currents in the DQD:s as a func- 
tion of the bias voltage Vi across DQDI. The current 
Ii through DQDi is determined by the populations of 
the eigenstates of the system and the effective tunnel- 
ing rates between the DQD and its right lead. The cur- 
rents Ii{Vi) and l2{Vi) are plotted for symmetric bias, 
f^Li = -M-Ri = eVi/2, and 0i = 62 = 7r/3 in Fig. g] 
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The key feature of both current-voltage characteristics is 
a thermally broadened onset at eVi = hu. In DQDl the 
onset occurs when the energies at which the electron can 
tunnel into and out of the DQD enters the bias window. 
The electron tunneling through DQDl will further excite 
the system to states where tunneling out of DQD2 be- 
comes possible. Hence the onset of the current I2 occurs 
at the same bias voltage Vi. In the limits ri/r2 ^ 1 
and ri/r2 S> 1 the solution to the ME, and hence the 
currents, can be obtained analytically (See Appendix C). 
Focusing on the non- locally induced current /2(^i)i for 
symmetric bias and eVi 3> ksT we obtain 



0.03 



i2 = ?^, ri/r2«i 



1-^7 



h = 



er27' 



[1 + 2 cot2(6'2)](l - 27 + 272) + r 



(4) 

Ti/Fa » 1, 



where 7 = fLi{fioj/2) cos'^{9i). From these expressions it 
is clear that the magnitude of the induced current can 
be made ^ eFi and ~ er2 in the limits ri/r2 <C 1 and 
ri/r2 ^ 1, respectively. From the plot of the high-bias 
current against the asymmetry factor ri/r2 for 9i = 62 
in Fig. 2] and from further investigation for 9i ^62, we 
find that the non-local effect is maximal for Fi ~ r2. We 
can can thus conclude that the non-locally induced cur- 
rent, qualitatively behaving as I2 ^ erir2/(ri + T2), is 
considerable for a large range of the parameters Fi, r2, 6*1 
and 02- 

To estimate the magnitude of the non-locally induced 
current we first recall that strong coupling implies the 
limit Fi <^ go for the tunneling rates. In recent single 
DQD-cavity experiments [13, fundamental frequen- 
cies !jj/2tt ^ lOGHz and DQD-cavity coupling strengths 
go ~ 500MHz were reported. This means that strong 
coupling requires F; < lOOMHz giving non-locally in- 
duced currents of the order I2 ^ O.lpA. Importantly, 
currents of this magnitude have been measured in DQD- 
cavity systems fSTj . 

To further put the magnitude of the non-locally in- 
duced current in perspective we briefly discuss non-local 
transport properties for the system with the DQD:s re- 
placed by single-level, or metallic dots. As follows from 
our Ref. [32[, the ME:s describing the evolution of 
these systems are explicitly dependent on the parame- 
ter Zq I Rq , typically much smaller than unity for regular 
transmission lines. Importantly, the ME:s show that the 
efi'ectivc tunneling rate into an empty dot or out of an 
occupied dot is independent on the system state to first 
order in Zq/Rq. It thus follows that the non-locally in- 
duced current will be proportional to {Zo/RqY, to first 
non- vanishing order. In a single-level dot- or metallic dot- 
cavity system the non-locally induced current will thus, 
in contrast to the current in a DQD-cavity system, be 
significantly suppressed. 
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FIG. 4: Upper panel: The non-locally induced current I2 as 
a function of eVi/hw for 61 = 62 — 7r/3, ksT — 0.05hw and 
different asymmetry factors ri/r2. The dashed black(green) 
curve shows the current for ri/r2 — 1/10(10) obtained from 
the analytical expression for Fi <^ r2(r2 ^ Fi) in Eq. ((5|. 
Inset shows Ii as a function of eVi/hw for the same asymme- 
try factors and mixing angles. Lower panel: The current I2 
above the onset voltage against the asymmetry factor ri/r2 
for different 6 = 61^62. 



B. Current correlations 



Having established that photon exchange between 
transport electrons in the spatially separated DQDis can 
result in a non-local current we, as the next natural step, 
investigate the mechanism behind this exchange. To this 
aim we consider the low-frequency correlations Sij be- 
tween the currents in DQDi and DQDj. Current cor- 
relations are known to provide information about e.g. 
the effective charge, interactions and statistical peorper- 
ties of the charge carriers (s^l- The correlations Sij can 
formally be obtained from a number-resolved version of 
the ME (See Appendix D). Focusing first on the cross- 
correlations 5*12 we plot in Fig. O the cross-correlation 
Fano- factor F12 = S12/ (eVliJ^) against the bias voltage 
Vi for different asymmetry factors. Similar to the nor- 
malized currents Ii and I2 the cross-correlations have an 
onset at eVi ~ hio. However, in contrast to the current 
I2 (but similar to Ii ) the cross-correlations have a strong 
dependence on the asymmetry factor F1/F2. This can be 
seen by considering the limits F1/F2 ^ 1 and F1/F2 ^ 1 
where analytical expressions can be obtained. Above the 
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onset voltage we get 



cos(gi)[l + cos4(gi)] 

-F12 = o(v/rvf^), Ts/ri « 1. 



(5) 



From these expressions we see that the currents in DQDl 
and DQD2 are manifestly positively correlated, F12 > 
(cos(0i) > 0), for T1/T2 ^ 1. The correlations are also 
strong, F12 ~ 1. In contrast, for ri/T2 ^ 1, the currents 
are essentially uncorrelatcd. The cross-over between the 
two regimes is shown in Fig. [5] for 61 = ^2- The strong, 
positive correlations appearing for Fi < F2 clearly show 
that tunneling through DQD2 is triggered by tunneling 
through DQDl. The qualitatively different system be- 
havior in the limits F1/F2 ^ 1 and F1/F2 S> 1, re- 
spectively is also manifested in the auto-correlations S'22- 
In Fig. [5] we see that the auto-correlation Fano-factor 
F22 = S22/{el2) above onset goes from a sub-Poissonian 
value, F22 < 1, to a super- Poissonian value, F22 > 1, as 
the asymmetry factor r2/ri is decreased from infinity to 
zero. This describes a transition from anti-bunching to 
bunching behavior of the transport electrons [s^l • 

To connect these findings to the properties of the pho- 
ton exchange we perform a careful investigation of the 
ME in the limit Fi <C F2. The processes contributing to 
to transport quantities to leading order in the asymme- 
try parameter F1/F2 are depicted in Fig. |6l From this 
scheme it is apparent that the states \D^0) and 152 0) 
will have occupations 0(1), while the other states have 
occupation C'(Fi/F2). The system will thus spend most 
of its time in the states \D^0) and I^^O) and will occa- 
sionally be excited out of this subspace by a tunneling 
event in DQDl, from the state I^^O) to any of the states 
{|D"1)}. The system can from here go back and forth 
between {IZJ^I)} and 15*20) an arbitrary number of times 
before relaxing to \D^0). This will occur on a time-scale 
~ I/F2. Each tunneling event in DQDl which excites the 
system to {|Z)"1)} will hence be followed by one or more 
tunneling events in DQD2 during a short time-window 
~ l/r2. The electrons in DQD2 are thus transported in 
cascades induced by randomly occurring tunneling events 
in DQDl with separation l/Fi ^ l/r2. This mecha- 
nisra, commonly referred to as dynamical channel block- 
ade [lll,!!^, explains both the positive cross-correlations, 
5i2 > 0, and the bunching of electrons in DQD2, F22 > 1. 
Importantly, each cascade in DQD2 is initiated by the 
emission of a single photon emitted by an electron tun- 
neling through DQDl. Our investigation thus supports 
the physical picture where single photons are exchanged 
between the DQD:s for Fi < r2. 
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FIG. 5: Upper panel: Cross-correlation Fano factor F12 = 
S12/ [ey/ I1I2) a function of eVi for different asymmetry fac- 
tors ri/r2 for mixing angles Oi — 62 — n/3 and tempera- 
ture fcsT = 0.05&J. Middle panel: Cross-correlation Fano 
factor above onset as a function of asymmetry parameter 
ri/r2 for 61—62—0 and /sbT <C Lower panel: Auto- 
correlation Fano factors F22 = S22/I2 above onset as a func- 
tion of asymmetry parameter for mixing angles Oi = O2 = 9 
and fcflT hu). 



C. Full counting statistics 



To obtain a complete picture of the elementary pro- 
cesses of the charge transport for Fi ^ F2 we consider 
the full transport statistics. The statistics is most clearly 



visualized via the cumulant generating function (CGF) 
F which can be obtained analytically above onset (See 
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a) 

tZ?"l) \D+l) 




|D°0) 

b) 

DQD2 

(2) (1) 




FIG. 6: (a) Schematic of processes in ME dP/dt = MP con- 
tributing to tranbsport quantities to first order in the asym- 
metry parameter Ti/r2- The tunnehng processes in DQDl 
(slow processes) are marked by blue arrows and the tunnel- 
ing processes in DQD2 (fast processes) are marked by red 
arrows, (b) Tunneling processes in DQD2. Process (1) and 
(2) describe the tunneling into and out of the DQD as the 
system goes back-and forth between ISiO) and {|Z)"1)}. Pro- 
cesses (3) and (4) describe the tunneling processes where the 
system relaxed from ISjO) to \D°0). 

Appendix D) and is given by 

-^(xi,X2) = -^|^l + cos2(0i) 

- ^sin''(6li) -f 4 cos2(6li)e*xi (sin2(6li) -f 0082(6*1 

y = (cos2(^?2) + sin2(02)e'x^) E„=o C^"'"'- (6) 

Here z = cos2(6'2)/[l + 0082(6*2)], zq ^ 1 + 0082(6*2) and 
Xi is the counting field for charge transfer in DQDi. 
To interpret the CGF we first consider the case when 
the charge transfer through DQD2 is not monitored, 
i.e. X2 — 0. Then y ~ 1 and the CGF reduces to 
the well known result for a single level dot with 
tunnehng rates Ti and Ti 0082(6*1) into and out of the 
dot, respectively. For X2 7^ 0, it is clear from the term 
sin^(0i) -|-oos^(6*i)j/ that tunneling events into DQDl are 
of two kinds. First, for events corresponding to the tran- 
sition IS'2) — )■ |-D°0), occurring with a rate oc sin^(6*i), 
there is no tunneling in DQD2. The second type of 



events, corresponding to transitions IS'20) — )• {|D"1)} 
and occurring with a rate oc cos^{9i), trigger tunneling 
in DQD2. This tunneling takes place as the system goes 
back to \D'^0) and is described by the function y = y{x2)- 
The different terms in y are given by the probabil- 
ities for all possible processes taking the system from 
to |D°0), weighted with counting field factors 
describing the respective charge transfer through DQD2. 
Common for all processes is that they start with the tran- 
sition {IL*"!)} ^ |S'?0) and end with jS'S'O) ^ |D°0). 
These two transitions give rise to the prefactor cos2(6*2) -I- 
sin^(6*2)e*>^2 = e*>^4cos2(6*2)e-*x=+sin2(6*2)], with the e*'^^ 
from the starting transition and cos2(6*2)e~*^2 -(- sin^(02) 
from the ending transition. As shown in Fig. 6, the end- 
ing transition can occur via tunneling from the left lead 
of DQD2 (probability sin^(0i), no counting field factor) 
or back from the right lead DQD2 (probability cos2(6*i), 
counting field factor e~''<^^). The sum term in y, running 
from zero to infinity, describes all possible back-and-forth 
transition \S'iQ) {\D°'l)] ^ |S'?0) the system can per- 
form, between the starting and ending transitions (see 
Fig 6). The n:th term in the sum thus corresponds to n 
transitions and hence n electrons being transferred across 
DQD2. This interpretation comes naturally when noting 
that 2;"/zo is the probability to return to the state \D"1) 
n times. The structure of the CGF, and in particular 
y, clearly shows that that electrons in DQD2 are trans- 
ported in cascades and that these cascades are triggered 
by single photons emitted by electrons tunneling through 
DQDl. 



IV. SPECTRAL FINE STRUCTURE AND 
ENTANGLEMENT 

It is interesting to investigate what qualitatively new 
physical effects come into play in the regime where the 
thermal broadening in the leads is much smaller than the 
DQD-cavity coupling strength, i.e. fc^T ^ hg^. In this 
regime the ME-description used in the previous section 
is no longer valid and we need to consider the full QME 
ofEq. dsn. 



A. Transport properties 

We first demonstrate that the structure on the scale 
^ in the spectrum of the generalized TC Hamilto- 
nian appear in the transport properties of the system 
for ksT < hgo. This fine structure, a manifestation of 
coherent electron-photon interaction, appear already in 
the average current and we therefore focus on this quan- 
tity. Moreover, we consider the simplest possible parame- 
ter regime by taking DQD:s with identical mixing angles 
6*1 = 02 = 0, lead-DQD tunnehng rates Fi = r2 = F 
and with left (right) leads having the same chemical po- 
tential, i.e. fi,^i = /i,y2 = fJ-u- This will give the same 
current / = /i = /a in DQDl and DQD2. This cur- 
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rent / is readily obtained from a numerical solution of 
Eq. (O (See Appendix D for details). In Fig. [7] we plot 
I{V) for symmetric bias voltage = —f^R ~ eV/2 and 
9 = for temperatures ranging from ksT ^ hgo down 
to ksT <^ hgo- We clearly see how the single step onset 
at eV = huj is split up into several smaller steps, spaced 
^ hgo as the temperature is decreased. These steps can 
directly be attributed to the structure of the spectrum of 
the generalized TC Hamiltonian. They are a consequence 
of eigenstates with energy splittings ~ hgo becoming pop- 
ulated at different bias voltages. It is here interesting to 
note that signatures of the JC spectrum was found in 
the frequency-dependent current auto-correlations in the 
transport through a system with only one DQD coupled 
to the cavity mode |24| . 



0.225 



l-H 




0.175- 



1.2 
eV/hro 



FIG. 7: The current I = h = h through the DQD:s as a 
function of voltage above onset for 6 = tt/3, go ~ O.lo; and 
different temperatures. 



the form in Eq. ([7]) the concurrence C{pr) takes on the 
simple form 



C = 



^{\PS - Pt\- 2^PePg,0} 



(8) 



To determine if entanglement can be induced by photon 
exchange we first consider the scheme in Fig. |8l dis- 
playing the lowest energy states with the TC-subspace 
well resolved. By noting that the lowest excited TC- 
states \D"'l) are written \D°l) = |50) and \D^l) ^ 
{l/^/2){\ggl) ± |T0)) it is clear that a selective popu- 
lation of any of the \D°'\) states would give an electronic 
state with a large singlet (5) or entangled triplet (T) 
component. To demonstrate such a selective population 
we choose bias voltages Vi and dot-level positions such 
that the chemical potentials jiiL and obey the rela- 
tions 

hhj f- huj 

— - V2hgo < iiLi < ^ - %o, 

hw hw , ,— , 

- — >Mffi > -^-(V2-l)fi5o. (9) 

For ksT <C hgo then only \D^1) of the excited TC-states 
becomes populated. 



1 000} 



|5J'0) 
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-■= 1 \D"l) 
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B. Transport-induced entanglement 

A natural question to ask when considering two cou- 
pled, spatially separated, DQD:s is to what extent their 
orbital degrees of freedom become entangled by the ex- 
change of cavity photons. The object of interest, describ- 
ing the properties of the electronic state with one electron 
in each DQD, is the reduced two-particle density ma- 
trix pr- The reduced density matrix, of dimension 4x4, 
is formally obtained by first projecting the total system 
density matrix p onto the TC subspace and then tracing 
out the photonic degrees of freedom. As follows from the 
structure of Eq. ([3|) and the TC-eigenstates (Appendix 
B), the reduced density matrix can be written as a sum of 
the four diagonal components in the singlet-triplet basis, 

Pr = Pa lag) {gg \+Pe\ ee) (ee \ + p s\S) {S\ + pt\T) (T | (7) 

where |5(T)) = {\eg) - {+) \ge))/^. The entanglement 
of Pr is conveniently quantified via the concurrence jssj , 
ranging from 1 for a maximally entangled state to for a 
non-entangled, separable state. For a density matrix on 



FIG. 8: Scheme of the lowest energy levels in the generalized 
TC Hamiltonian. The red arrows show the active transitions 
for fcsT ^ hgo and for a bias configuration such that the 
chemical potentials satisfy the conditions in Eq ((9]). 



For the chosen parameters, as seen in Fig. [8l only five 
states of the generalized TC-model contributes to trans- 
port. For this case the QME can be solved exactly. Im- 
portantly, the steady-state solution gives a reduced two- 
particle density matrix with only pg and px nonzero. 
As is clear from Eq. ([HJ the resulting concurrence is finite. 
For a symmetric parameter setting, i.e. for 9i ^ 62 = 0, 
Fi = F2 = F and fi^i = /i^2 = Mi^ this concurrence is 
given by 



C = 



2[cos4(6') -|-sin'^(6l)] 



(10) 



showing that the concurrence can reach up to C = 1/2 
for 9 <^ 1. We stress that for the chosen parameters there 
is no bound on 9 in order to have a well defined solution 
to the QME. 
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Having confirmed the existence of large entanglement 
C < 1/2 we consider the effect of finite temperature and 
modified bias voltage V = Vi ^ ¥2- Solving numerically 
the QME, the resulting concurrence C{V) is plotted in 
Fig. [9] for different temperatures. For low temperatures 
ksT <^ hgQ the entanglement has an onset when \D^1) 
is populated, with the concurrence given by Eq. (|10p . 
Increasing the voltage further \D^1) is populated as well, 
decreasing the concurrence due to the finite probability 
for both entangled triplet, pt and singlet, ps, electronic 
states, clear from Eq. For even larger bias all TC- 
states have finite population and the entanglement 

disappears. From Fig. |9] it is also clear that increasing 
the temperature smears the C{V) curve and successively 
suppresses the entanglement, reaching a separable state 
at ksT hgo/2. 
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FIG. 9: Concurrence C as a function of bias voltage for dif- 
ferent temperatures for go ~ O.Ilj, pl + ^j? = — 0.12fitt; and 
91=92 = n/6. 



We can thus conclude that both in the high-bias and 
high-temperature regimes, where the system can be de- 
scribed by a ME, the entanglement is zero. To clarify the 
generality of this observation we investigated the concur- 
rence in all regimes where we could solve the ME an- 
alytically (see Appendix C). In these regimes we could 
formally prove the absence of entanglement. Moreover 
we considered the concurrence obtained numerically for 
a broad range of other system parameters in the ME- 
regime but did not find any entanglement. We thus 
conclude that it is highly probable that electrons in 
the two DQD:s can only be entangled for temperatures 
ksT <^ hgo, in biasing regimes where TC-states with the 
same number of excitations are selectively populated. 

We stress that our investigation mainly aims at demon- 
strating the existence of entanglement. We do not ana- 
lyze how or even if it can be detected by transport mea- 
surements. Moreover, we have not made a full investi- 
gation over the entire parameter space to identify the 
regime with the largest entanglement. 



V. DEPHASING AND RELAXATION EFFECTS 

So far we have neglected dephasing and relaxation ef- 
fects in the DQD:s as well as loss of cavity photons. From 
the recent single DQD experiments [13, [1^ it is clear that 
the dephasing rate F^i is much larger than the rates Tn 
and K for relaxation and cavity loss, respectively. We 
thus focus on the effect of dephasing on the results pre- 
sented above. 

Dephasing can qualitatively be accounted for by 
adding a term [s^ 



Cd[p] 



2 



[2L\pU-L\Up-pL\l^ (11 



1=1,2 



(f d 



to the Liouvillian in Eq. ([3]). Here Li = d\idei — d'g 
and the dephasing is taken independent, with the same 
rate Td, for the two DQD:s. An investigation of de- 
phasing in all parameter regimes is beyond the scope of 
the present article. However, we stress that for strong 
dephasing. Yd ^ go, coherent superpositions between 
excited and ground states in the DQD:s are suppressed. 
The steady state solution of Eqs. (O and pT|) is diago- 
nal in the basis of the DQD-cavity product states \^i£,2p), 
with = 0, g,e and p the number of photons. As a con- 
sequence, electrons and photons are decoupled and the 
non-local transport effects as well as the DQD entangle- 
ment appearing in the regime hgo 3> ksT are suppressed. 

Importantly, in both single DQD-experiments [13, 
the dephasing is found to be strong, with Tjj ^ 1 GHz, 
substantially larger than the coupling strength 5o/27r ^ 
100 MHz. It is thus necessary to consider ways to increase 
go and/or suppress Td, in order to approach the strong 
coupling limit go ^ To where the non-local effects dis- 
cussed above are fully developed. First and foremost, the 
coupling go can be increased substantially by increasing 
the fundamental frequency {go oc cu), simply by making a 
shorter cavity. Importantly, since we consider an isolated 
cavity, a; is not limited by requirements of an external mi- 
crowave circuitry. The limit is instead set by the energy 
gap of the superconducting cavity material, of the or- 
der of hundreds of GHz for large gap superconductors as 
e.g. Nb (w/27r « 10 GHz in [33, [S^]). Second, unconven- 
tional transmission line cavities, with a central conductor 
consisting of e.g. Josephson junctions or SQUIDS, (60l - 
[5^ can have characteristic impedances Zo ^ lkf2. This 
gives coupling strengths go ~ O.lw, one order of mag- 
nitude larger than for conventional transmission lines. 
Third, since Refs. [13, [1^ are the first experiments on 
DQD:s in cavities, there is probably room for for optimiz- 
ing the circuit design, further suppressing the dephasing. 
Taken together, the strong coupling limit of our proposal 
is arguably within reach experimentally. Moreover, the 
relaxation rate in [1^ was estimated to F^ 100 MHz, 
one order of magnitude smaller than the dephasing rate. 
Hence, in the strong coupling limit go S> Fo, relaxation 
is expected to be negligible. In addition, the cavity loss 
rate k in (stI. [39| was already much smaller than go, sug- 
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gesting that cavity loss can safely be neglected. 



VI. CONCLUSION 

In conclusion we have theoretically investigated the 
non-local transport properties of a DQD-cavity system. 
We have found that the photons emitted by electrons 
tunneling in one DQD can assist transport of electrons 
through the other DQD. giving a strong non-locally in- 
duced current and large cross-correlations between cur- 
rents in the two DQD:s. Moreover, in the low temper- 
ature regime, fc^T <^ hgo, we have demonstrated that 
signatures of the TC-spectrum will appear in the I-V 
characterstics and that the orbital degrees of freedom of 
electrons in the two DQD:s can become entangled. Im- 
portantly, our work provide a theoretical framework for 
investigations of non-local electronic transport proper- 
ties in cavity -coupled nanoscale conductors. The anal- 
ysis can readily be modified to study transport through 
other nanoscopic two-level systems coupled to cavities, 
e.g. superconducting single electron transistors [13, H^] 
and spin qubits [2l|, [sl, IH, [63| 
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Appendix A 

We here derive the Hamiltonian for the DQD- 
transmission line cavity system. The first step is to 
describe the DQD:s within the standard [13] constant- 
interaction model. Then only the excess dot charges will 
interact capacitively with the cavity. This means that the 
total Hamiltonian of the system, Hs, becomes the sum 
of the Hamiltonian for the orbital states of the DQD:s, 
Ho, and the Hamiltonian for the cavity, the dot charges 
and their interactions. He- For the orbital part Hq we 
consider DQD:s formed by two tunnel coupled quantum 
dots with a single active spin-degenerate level in each dot. 
The orbital part of the Hamiltonian, in the localized basis 
of the DQD:s, then has the form 



dmdLi) 
(12) 



with Ai being the energy difference between the bare en- 
ergies of the orbitals in the left and right dot of DQDi. 
We recall that tLm denotes the interdot tunneling am- 
plitude of DQDi and note that the creations operators 
d^Lvd^Ri related to the eigenbasis creation operators 
creation operators according to 



dh- 



sv[i{9i)d 



cos 



= cos{0i)dl^ + sin 



mdl 

mdi- 



(13) 



The Hamiltonian He is derived within the frame- 
work of circuit QED. Following the procedure of Refs. 
[3^ [gsI |66| we start from the classical Lagrangian of a 
circuit representation of the system, including the capac- 
itances of the dots. The transmission line is modelled 
by a single LC-circuit. This will describe the physics 
of one finite-frequency (the fundamental mode) |67| and 
the zero-frequency mode. The circuit diagram is shown 
in Fig. [ini where Cci, Cg^i, Clri and Ci, Vg^i denotes 
the capcitances and gate voltages of DQDi and Lq, Cq 
are the total inductance and total capacitance to ground 
of the central conductor. 
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FIG. 10: Diagram of the circuit describing the transmission 
line cavity, the dot charges and their interactions. In DQDI 
(DQD2) the right (left) dot is coupled capacitively to the 
central conductor of the transmission line, modelled by a LC- 
circuit and to two gate electrodes. The nodes vi correspond 
to the dots while node 1 and 2 correspond to the endpoints 
of the transmission line. 

The Lagrangian of the circuit is given by 



+ 



CLRi{(t>Li — <t)R.if , X^jy Cgvi{4>„i — Vg^i)'^ 



2L, 







(14) 



where (f>^i and </),; denote the phases of node vi and 
i, respectively (see Fig. [TU]). The zero- and finite- 
frequency normal modes, describing the electrostatics 
and -dynamics of the circuit, respectively, are obtained 
from the Euler-Lagrange equations. By rewriting the La- 
grangian in terms of these normal modes, performing a 
Legendre transformation and a canonical quantization of 
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the fundamental mode a quantum Hamiltonian, He, is 
obtained. By further writing the exeess charges of the 
dots as e{dl^d,yi — Ug^i), where n^^i denotes the gate- 
induced charges, we get He = Hj + Hoia + Hjjii, with 

Hj ~ hijja^a 

e^{dld _-ng^,)^ ^ ^^^^.^ ^ ^^^^^j^^ _ ^^^^^ 



E 



Hdio = E Ui^it^^idiidi^i - ngi,i){d^^2df,2 - "-9^2) 



tical couphng strengths gi = 52 = go, considered in the 
main text, the mixing angles 61 and 62 cannot be tuned 
independently for fixed Xli and Xm. 



Appendix B 

We here give the explicit form for the eigcnstates with 
a finite number of excitations in the TC-subspace and all 
eigenencrgies of Hs ■ The former are given by 



Hon - 



\D^1) = \S0), |,,±i)= 1-951) ±1-^0) 



V2 



(16) 



(15) 



and for n > 2 



In the total Hamiltonian Hs = He + Hq, the part 
Ho + Hj has the standard form for a few-level dot 
system linearly coupled to a bosonic mode (see e.g. 
Ref . [68t ] . with C^i and A^/i denoting the effective self- 
capacitances of the dots and the coupling strengths be- 
tween the photons of the fundamental mode and the 
dot charges, respectively. The part Hdio contains cross 
terms, ^ d]^j^d^idj^2d-tj.2, describing direct non-local cou- 
pling between the dot charges. The coupling strengths, 
Uvifi2, depend parametrically only on the capacitances 
of the circuit in Fig. [TUl i.e. not on the inductance Lq, 
and would thus remain unchanged if this inductance was 
short-circuited. This means that Hdiq describes purely 
electrostatic, or capacitive, coupling. The part Hun, 
just as Hdio, describes direct non-local coupling between 
the dot charges. However, in contrast to Hdio the cou- 
pling strengths in Hdii, X^iX^a/fu^, depend parametri- 
cally on the inductance Lq. This part is therefore elec- 
trodynamic. 

As the next step we motivate the approximations lead- 
ing from Eqs. ((T^ and (fT5|) to Eq. (P), under the 
conditions described in the main text. To do this we 



\D+-n) 



\/n — 1 \ggn) — ^Jn \een — 2) 



, \ggn) + \/n — 1 \een - 2) \Tn — 1) 



v/2(2n-l) 



± 



\D-+n) = \Sn-l) , 



V2 ' 
(17) 



with \S{T)n) = {\egn) - {+) \gen))/^/2. The eigenener- 
gies are 

eoon = fi^(n- + 1) (18) 
for the HO subspace, 

£500 = fiw/2, e5±„ = hio{n + 1/2) ± ^/nhgo, n > 1 



(19) 



for the JC subspaces, and 

e_D"o = 0: = eD±i = hio ± \/2hga. 



(20) 



use the re lation s A^; cx yJZo/Rqhhj {Zq = y/Lo/Co), for the TC subspace. 
w ~ 1/^/LqCq and U^i^2 <x. e^/Co for the parame- 
ters in He- First, noting that Zq <C Rq and hence 

Xyijhbj <C 1 justifies a rotating- wave approximation. Appendix C 

which amounts to neglecting all terms of OiXvijuS) in 
Hs, e.g. the counter-rotating terms. Second, for the 
DQD:s resonant with the cavity mode the direct capac- 
itive interaction will scale as J7;yi^2/Ay(p)i ~ \/ Zq/ Rq 
and can thus be neglected. Moreover, had we consid- 
ered a full description for the transmission-line DQD cir- 
cuit, including all of the cavity modes, the higher fre- 
quency modes would be off-resonant with a detuning 
/S.E > fko. The corrections due to this off-resonant in- 
teraction would then scale as Xyi/tnj and therefore be 
negligible. It should be noted that for the resonance con- 
dition to hold the tunneling amplitudes and detunings 
between the left and right dot-o rbitals the DQ D:s must 
be chosen such that |Aj| = 2yJ {hwY — iLm- 1™" 
ther point out that the DQDi-cavity coupling strengths 
is gi = sm-{29i){Xiii — XLi)/2. Thus, for the case of iden- 



In this Appendix we explain how the QME in Eq. ([3]) 
can be reduced to a ME in the limit hg^ ^ kgT, give 
the explicit form of the ME and solve it in three lim- 
iting cases. We start by pointing out that in the sec- 
ular regime ^ g^, considered here, only coherences 

between degenerate states, i.e. | H n) and | — -I-71) 

need to be taken into account in the QME. Moreover, 
only the diagonal elements {Sfn\p\S°'n) couple to the 
coherences. As pointed out in the text, for hg^ <^ ksT, 
the QME becomes independent on g^. This introduces 
additional symmetries in the QME, with two important 
consequences: (1) The coherences (-1 — n\p\ — \-n) and 
( — — n) couple with opposite signs to {Sfn\p\Sfn) 
and {S~ n\p\S~ n) . (2) For several pairs of diagonal ele- 
ments of p, only the sums of the elements couple to the 
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other diagonal elements. In particular, this holds for the 
sum {S^n\p\S^n) + {S~n\p\S~n), to which the coher- 
ences, according to (1), do not contribute. As a result 
of (1) and (2), the coherences decouple from the diago- 
nal elements of the QME, allowing us to reduce it to a 
standard ME. 

To write the explicit form of the ME it is convenient 
to first introduce a shorthand notation for the diagonal 
elements of p, i.e. the probabilities for the eigenstates of 
Hs- The probabilities, or the sums of probabilities, for 
states the HO, JC and TC subspaces are denoted by 

Po"o - (00n|p|00n) , 

Pg = <5„o {S^0\P\S^0) + (1 - Sno) J2 iS>\p\Srn) 

a = ± 

(21) 

and 

= ,5„o (i?°0|p|i?°0) + 5„i ^ {D"l\p\D"l} 

+ (1 - (S„o - Snl) J2 (^""^l/5|^""^) , 

Pho = {DH\p\D"l) , 

P^+_(_+, = {D+-(-+'>n\p\D+-^-+'>n) , n>2, (22) 



respectively. By further introducing vectors Px = 
(P^ P^ P|.....)^, with X = 00,51,5*2, containing the 
probabilities for states with one or both DQD:s unoccu- 
pied and vectors P_d = (P|, Pjj....)^ and Pjj+-(-+) = 
(P^^_(_^) P^^^_^_^y...)'^ containing the probabilities 
for states with both DQD:s occupied, the ME can be 
written 
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(23) 



The submatrices in M above and below the diagonal are here given by 



J=0,1 j 



jj _ 4n 3 - 2j 



n + j 



4(2rH-l) ' 2n^V 
(M|;o),„ ^ ^i„^o^ (M|^°i)i„ = (5i™gi/2, (M|^i)i„ = <5i™&^/4(A/^-o)„i = (5„i^o_ 



1 

4 



(24) 



and 



(M00)„„ = ^ 



j=04 



2 — <^nl'5il 



D 



(4n + 3 + 2j) 

Vnj 



J=0,1 

n+j _ ,,_D-+ ^ 1 
2' 



"J 4(2n + l + 2j)' '""^ 2n-hl-h2j' 
(M^'oi)«i = <5ni^r/2 + <5„2^- /2, (M^iU = SniQV^ + <5«23e,V4, 



(25) 
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respectively. Here = f ,g J,,(-fi£^/2), G} = j:^f,,,U,{nLj/2), Gf = f ,gJ,,(-Sw/2), G} = 
X^i/ -ri^ei/i/i(^w/2) and we use the i-index convention 1 = 2 and 2 = 1. The submatrices on the diagonal in M 
are diagonal with elements such that the sum of every column in M is zero. This structure ensures that the ME 
conserves probability. 



We find the steady-state solution to the ME in Eq. 
analytically in three limiting cases: 

(i) For 6*1 = 6*2, symmetric bias voltages eVi,eV2 ^ 
applied across both DQDl and DQD2. 

(m) For the bias condition of Sec. IIIII with Fi 3> F2 and 

symmetric bias voltage across DQDl. Here we calculate 

the distribution to zeroth order in F2/ri. 

{iii) For the bias condition of Sec. IIIII with F2 ^ Fi and 

symmetric bias voltage across DQDl. The distribution 

is here calculated to first order in F1/F2. 

The solution for case (i) is used to derive the stability 

condition in Sec. IIIBI while the solutions for the cases 

(ii) and (iii) are used to obtain the analytical expression 
for the current in DQD2, correct to first order in F2/F1 
and F1/F2, respectively. (See SeclIIip. 

To find the solution of the ME in case (i) we use a gen- 
eral property of the ME. This property states that Pg^ 
and Pg^ will couple only to the probabilities for states 
with both DQD:s unoccupied having n or n -|- 1 excita- 
tions and for states with both DQD:s occupied having 
n or n — 1 excitations (See Fig. [2|). In turn the prob- 
abilities for these states wiU couple to Pg^'' and Pg^'" , 
with k = —1,0, 1. Two coupled second order difference 
equations are thus obtained for the probabilities in Ps-^ 
and Ps2- For the special conditions (i) these difference 
equations become particularly simply. Introducing the 
vector P„ = {Pg^ Ps^)'^ these equation be written 



sin''(6l)MoPi 



sin''(6l)MiP2 

- 2cos''(6l)AfoPo, 
sin4(0)A/„P„+i = 



2cos'*(6')A/o + [1 
cos (6')Afi -f cos^(6')A/o 



cos^{e)]B) Po 



Pi 



cos4(^?)M„ + sin4(0)A/„_i + ^ 



-cos4(6')M„_iP„_i, n>2 
with 



(26) 



B = F1F2 



1 -1 
-1 1 



and 



rl/2 + Tlbn FiF2(l/2 + 6„ 
FiF2(l/2 + 6„) Fi/2 + F2&„ 

(4n-t-3)(4n+l) 
' 4(2n-|-l)[4n+l + 2cos2(0)] 
n{n + 1) 



(27) 



2(2n-t- l)[n-t-sin^ (61)] 



Together with the condition P„ — > for 



(l26l) has the solution P,- 



cot 



4n+2( 



(28) 



00 Eq. 



the probabilities for the states with both DQD:s occu- 
pied and unoccupied, respectively. Requiring that the 
solution to Eq. (1^51) is normalized we get 



P^l = cot4(«+i)(0)a;o, Pi - cot2(0)xo, 
P^^ = 2cot^"+2(6i)a;o, n>l 
P^=XQ, Ph = 2PAo = 2cot4(^).To, 
Pi5 = 2P^+_(_+, = 2 cot4"(0)xo, n>2, 



(29) 



witha-o = [l-cot**(6')]/[l + cot2(6')+cot'^(6')]^ It is clear 
from Eq. ()29p that a well defined solution, or equiva- 
lently a solution with non-infinite mean number of ex- 
citations, exists if > 7r/4. On physical grounds we 
argue that there exists a more general stability condition 
applying also for mixing angles 9i ^ 62. We start by 
considering the case when 0i = 6*2 > 7r/4 and one of the 
mixing angles is increased. From the discussion below 
Eq. ([3]) it is clear that this will increase absorption rela- 
tive to emission of photons by tunneling electrons. The 
mean number of excitations will then be decreased and 
the distribution must therfore still be convergent. Since 
all pairs of mixing angles ^1,^2 > 7r/4 can be reached this 
way it follows that a well defined solution exists for all of 
them. The mean number of excitations will also decrease 
if the bias voltage is decreased. The stability condition, 
6*1, 02 > ''i'/4, must therefore also hold for finite bias volt- 
ages. This conclusion is further supported by numerical 
investigations. 

In the limiting case {ii) the relation between proba- 
bilities for states with DQD2 unoccupied, e.g. P^^ and 
Pg\, are entirely determined by the tunneling in DQDl. 
Similarly he tunneling in DQDl entirely determines the 
the relation between probabilities for states with DQD2 
occupied. The tunneling in DQD2 only effects the total 
probability for DQD2 being occupied. The main steps 
in the solution of the ME in this limit are most clearly 
visualized by rewriting MP = as 



M(°)+F2/riM, 



(1) 
00 



-"^11 



F2/riAliW \ / P20 

F2/riAfiW i I P21 



M 



(30) 



where the vectors P20 and P21 contain the probabili- 
ties for states with DQD2 unoccupied and occupied, re- 
spectively. The matrix M further contains the transition 
rates in Af divided by Fi. The starting point of the 
derivation is to note that for F2 = the matrix is block- 



{e){l l)'''Pi"i, inde- diagonal, with blocks Af^\ and that Aeii^Ml^') 



r(0)^ 



0. 



pendent on Fi and F2. This solution is then used to find It then follows that the eigenvalue zero of the matrix 
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M is doubly degenerate for r2 = 0. To find the so- 
lution P(°) to Eq. (1201) to zeroth order in Tz/ri, i.e. 
the limit of P as T2/T1 — > 0, wc must therefore apply 
degenerate perturbation theory generalized to ME ma- 
trices. The first step in this procedure is to find the 
two linearly independent solutions to Eq. (l30| . i.e. the 
eigenvectors corresponding to the eigenvalue zero of M. 
Setting r2 = in the equation the two linearly inde- 
pendent eigenvectors acquire the forms Pq = (P20'' 



MePe = 0. We get 



and Pi = (0 P^i^)^, where P^q^ and P^i^ fulfill the 
equations -A/qq''P2q^ = and M^^^P^^ = 0, respectively. 
Importantly these equations can readily be solved ana- 
lytically as they give difference equations similar to Eq. 
(P^l) . The next step is to express P'-^^ as a normalized 
linear combination of these vectors, i.e. cqPo + ciPi, to 
which P tend as r2/ri — ^ 0, . To do this we define the 
projectors 7^ Po(?7o 0)-hPi(0 Ui) and Q= l-V. 
Here {Un 0) and (0 Ui), with = (1 1 ...), are the 
left eigenvectors to the eigenvalue zero in M normalized 
so that UiP2i = 1. Then by applying the steps presented 
in Ref. [H the equation VM'^^^rP^^^ = for P(o) is 
obtained. In turn this equation gives 



(0) 



,(0) 



C0C/0M00P20 + C1U1M01P21 = 0, 



(31) 



which together with the normalization condition cq+ci = 
1. determines co and ci. The different element in P'*'^ 
can then be written 



^Jo = [7/(l-7)r+% 
P^\=(2-M[7/(l-7)r5o) 
P^, = (2-M[7/(l-7)r+'5i 



pO 



Cl, 



pi 



2P»|,=2[7/(l-7r5i 



pn 9 pn 



.(_+, =2[7/(l-7)]2"Si, n>2, (32) 



where 



Co = 



'77 



Cl 



(l- 27 + 272)[l + 2cot2(0i)] +72' 

7?(l-7^)[l + 2cot2(gi)] 
(l-7 + 72)[l + 2cot2(0i)]+72' 



(33) 



and 7] = (1 — 27) /(I — 7 + 72). With the probabilities in 
Eq. (132]) and the rates between DQD2 and its right lead 
for the corresponding states it is straightforward to derive 
the second expression for the current through DQD2 in 
Eq. 

To solve the ME in case (iii) we first recall, from the 
main text, that only the state and transitions depicted 
in Fig. 101 contribute to the transport quantities to first 
order in ri/r2. To find the probability distribution to 
P first order in ri/r2 only the additional transitions 
{|D"i)} — )■ 1 5*20) need to be taken into account. Thus, 
the distribution is obtained by solving an effective ME 



^ ^ _ ri 2[2 + cos2(0i)] 72 
1 + 7 r2 sin2(0i) (1 + 7)2 
Fi [l + cos2(gi)]72(l-7) 
r2 sin2(0i) (1+7)' 

7 Fi2[2 + cos2((?i)] 73 



p; 



^= 1 + 7 sin2(6li) (1 + 7)^ 

Fi [l + cos2(0i)] 7^(1-7) 



sin2(0i) (1+7)2 
Pd = 2Pi30 = 



1 _ 2Fi l + cos2(6'i) 72 



F2 sin2(0i) (1 + 7) 



r2(l+7)' 



(34) 



These probabilities are used to obtain the first expression 
for current through DQD2 in Eq. ©. 



Appendix D 

In this Appendix we present the derivation of cur- 
rents, current correlations as well as the full statistics 
of charge transfer accross the DQD:s. Following the 
procedure of Refs. [zO, [zll we rewrite the QME of 
Eq. (IHI in the n-resolved form and Fourier transform it 
with respect to the number of electrons having tunneled 
through DQDl and DQD2. The QME then transforms 
to dp/dt = £(xi,X2)/3, where the counting fields xi and 
X2 are the conjugate variables to the number of electrons 
having tunneled through DQDl and DQD2. The eigen- 
value of £(xi,X2) tending to zero as Xi>X2 — >■ is the 
long time limit cumulant generating function ^(xijX2)- 
The currents and the noise are obtained from the first and 
second derivatives of J^{xitX2), i-e. h = edi^-T\xi=x2=o 



and Sij = e'^di^-di^- J- 1^-^=^2=0 ■ These quantities can 
conveniently be accessed via the eigenvalue problem 
'C(xi,X2)[/5(xi,X2)] = ^{xi,X2)p{Xi,X2)- In the present 
paper this full QME-approach is used only to calculate 
the current in Fig. H 

In the ME limit the eigenvalue problem above reads 
^^(xi,X2)P(xi,X2) = -^(Xi,X2)P(xi,X2)- We use this 
equation to calculate the noise plotted in Fig. HI We 
also use the equation to obtain the CGF of Eq. i.e. 
the CGF to first order in Fi/r2 for the bias condition of 
DQD2 described in Sec. Ulllwith DQDl in the high-bias 
regime. To do this we write det[Af (xij X2) — -^(xi? X2)] = 
as 



det 



M,{xuX2)-J' M,^e(Xl,X2) 
Me^.(Xi,X2) Af.(xi,X2)--^ 



= 0. 



(35) 



Here Me-^r, Mr-^e and Mr are matrices describing the 
transitions to, from and within the subspace of states 
not included in the Fig. l6l Importantly, {Me^r)nm oc 
0(ri/r2) and {Mr^e)nni oc 0(1). This means that 



det[Me(xi,X2)--F(i)(xi,X2)] =0, 



(36) 
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where -^"^"'^■'(Xij X2) is the CGF to first order in ri/r2. for J'^^) then gives Eq. 
Dropping terms of C'[(ri/r2)^] in Eq. ([55)) and solving 
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